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ABSTRACT 

Previous work on the quasar mass-luminosity plane indicates the possibility that 
quasars of the same central black hole mass might follow a common evolutionary 
track, independent of the properties of the host galaxy. We consider two simple mod- 
els for the evolution of individual quasars. Requiring these tracks to lie within the 
observed quasar locus at all redshifts strongly constrains the model parameters, but 
does allow some solutions. These solutions include a family of tracks with similar shape 
but different initial masses that might match the observed quasar distributions at all 
redshifts z < 2.0. This family of solutions is characterized by short (1-2 Gyr) lifetimes, 
a duty cycle in which the quasar is on at least 25% of the time, and a rapid decline in 
Eddington ratio, perhaps with LEdd oc t~ 6 or steeper. 

Key words: black hole physics — galaxies: evolution — galaxies: nuclei — quasars: 
general — accretion, accretion discs 



1 INTRODUCTION 



Supermassive black holes (SMBH), are found at the cen- 
ters of nearly every galaxy where there have been sen- 
sitive searches. These SMBH are the end products of a 
long lifecycle comprising a seeding mechanism, a growth 
stage and, finally, the present state with nearly zero growth. 
' At various points during this cycle, individual SMBH are 
viewed in many states, including as Type 1 and Type 2 
quasa rs, Seyfert ga laxies, and quiescence. The Soltan argu- 
ment (|Soltanlll982l ) indicates that most of the final SMBH 
mass is gained through luminous accretion while the SMBH 
and its surrounding region are in a quasar phase. It has 
recently been demonstrated that quasar masses and lu- 
minosities are more strongly con strained at each redshift 
than e xpected from simple models (|Steinhardt fc Elvisll2O10l . 
l2009al ). In this work, we use these constraints as a probe of 
SMBH accretion tracks to produce limits on possible feed- 
back mechanisms. 

These constraints include new boundaries on quasar lu- 
minosities on both the bright and faint ends at 0.8 < z < 
4.1. Several examples of quasar distributions in the mass- 
luminosity plane are shown in §[31 and on the high-mass end, 
the entire quasar distribution lies above the Sloan Digital 
Sky Survey (SDSS) detection limit, bo unded on the bright 
end by a sub-Eddington boundary jSteinhardt fe Elvis! 
2010) and on the faint end by high -mass, low-luminosity 
boundary (|Steinhardt fc Elvis! l2009al l. At these redshifts, 



the entire range of quasar luminosities for a given mass 
spans less than 1 dex ever ywhere that it can be measured 
jSteinhardt fc Elvis! 1200931). Given the ~ 0.3 dex variabil- 
itv of quasars (|de Vries et al.l 120051 ). the intrinsic width of 
this distribution must be quite narrow. For fixed mass and 
redshift, the luminosity distribution is sharply peaked at 
one, characteristic luminosity. This luminosity drops signif- 
icantly below the Eddington limit for larger central black 
hole masses at each redshift. The lowest-mass quasars at 
each redshift can typically approach their Eddington lumi- 
nosity, while the highest-mass quasars at each redshift typ- 
ically are below ~ 10% of Eddington. At 0.2 < z < 0.8, the 
SDSS detection limit is within 1 dex of the bright-end, sub- 
Eddington boundary, so we cannot test whether the width 
of the entire luminosity range is less than 1 dex. However, 
there is again a sharp peak at one, characteristic luminosity 
and we might expect that with improved detection, quasars 
at 0.2 < z < 0.8 would exhibit identical properties to those 
at higher redshifts. The observed characteristic luminosity 
evolves with mass at fixed redshift and evolves with redshift 
at fixed mass. 

The remarkable implication of these constraints is that 
there is a characteristic luminosity for quasar accretion de- 
pendent only upon the central SMBH mass and redshift. 
In particular, even though we would expect a large sample 
of supermassive black holes of a given mass and redshift 
to lie within host galaxies of a range of different virializa- 
tion times, morphologies, star formation rates and merger 
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histories, the corresponding quasars have a common char- 
acteristic luminosity. It is also unexpected that the char- 
acteristic luminosity evolves with redshift but also remains 
narrow at fixed z. This combination appears to require some 
non-trivial evolution of accretion rates that occurs simulta- 
neously in different host galaxies. A discussion of physical 
models capable of producing this behavior is left for future 
work. In this paper, we use the evolution of these narrow lu- 
minosity ranges to constrain quasar accretion histories, i.e., 
tracks taken in the mass-luminosity plane. 

Recent work has investigated the possibility that 
virial mass estimates ma y cont ai n a mass-de pe ndent 
bias (lOnken fc Kollmeierl 120081 ; iRisaliti et all 120091 ; 
iRafiee fc Haul bold ). These proposed mass-dependent 
biases all have th e property that black holes measured in 
Sh en et al.l (|2008T ) to have similar masses also have similar 
masses after the corrections are applied (and black holes 
measured to have dissimilar masses have dissimilar masses 
after correction). As a result, the luminosity distribution at 
fixed mass and redshift continues to be equally narrow after 
correcting for any of these biases, and evolutionary tracks 
will yield comparably restrictive, though different, bounds 
to those described in this work. 

A bias that takes the true M — L distribution and maps 
it into one that fills less of the estimated M — L parameter 
space will affect the tracks that one derives (and perhaps 
the validity of the method for a strong enough bias). This 
method would produce less restrictive bounds under a cor- 
rection moving objects mea sured to ha v e sub stantially dif- 
ferent black hole masses in IShen et al.l (|2008l ) to the same 
corrected masses, and thus widening the luminosity distri- 
bution at fixed M and z. 

The quasars with lowest central black hole mass at each 
redshift can approach their Eddington luminosity. As they 
accrete, their mass increases, and as quasars move both to 
higher masses and later times, their characteristic Eddington 
ratio decreases. The hi gh-mass boundary is also s harp and 
evolving with redshift ( Steinhardt fc Elvisl l2009ah . so that 
most SMBH survive as quasars until they reach this bound- 
ary. The sharp, characteristic luminosity at fixed mass and 
redshift, and its seeming independence of host galaxy pa- 
rameters, imply that these surviving quasars take similar 
paths in mass and luminosity as a function of time from 
when they first appear near L E dd until turnoff. 

While for each individual quasar there is only one snap- 
shot available showing its mass and luminosity at one red- 
shift, we can reconstruct possible histories for that quasar 
by requiring that its evolution would, at each redshift, pro- 
duce a mass and luminosity lying within the observed quasar 
locus in the M — L plane, turning off when the quasar has 
become one of the highest-mass SMBH at some redshift. We 
find that the specific path taken is sensitive to the SMBH 
accretion history, and is sensitive to feedback between the 
SMBH and its environment. In § [3] we develop two simple, 
parametrized models f or SM BH accretion history following 
Hopk ms fc Hernquistl <|2009n and show that the resulting 
quasar tracks tightly constrain their input parameters. In 
§[4] we fit these models to the SDSS quasar mass- luminosity 
plane and consider the properties that seem to be required 
of these tracks. We discuss the implications of these results 
in §0 



Table 1. Locations of the sub-Eddington boundaries (SEB) and 
high-mass, low-luminosity boundaries (HMLLB) at 0.2 < z < 2.0, 
where M = log M BH /M Q . 



z 


SEB (logL max ) 


HMLLB (logL min ) 


H/3 


0.2-0.4 
0.4-0.6 
0.6-0.8 


0.37 M + 42.66 
0.45 M + 42.17 
0.61 M + 41.14 


detection-limited 
detection-limited 
detection-limited 


Mgll 



0.8-1.0 


0.67 M - 


t- 40.66 


0.54 M - 


- 40.97 


1.0-1.2 


0.67 M - 


t- 40.78 


0.35 M - 


- 42.77 


1.2-1.4 


0.73 M - 


h 40.30 


0.44 M - 


- 42.01 


1.4-1.6 


0.68 M - 


h 40.73 


0.47 M - 


- 41.80 


1.6-1.8 


0.51 M - 


t- 42.40 


0.44 M - 


- 42.20 


1.8-2.0 


0.42 M - 


h 43.24 


0.31 M - 


- 43.50 



2 MASS-LUMINOSITY PLANE BOUNDARIES 

As shown in Papers I (ISteinhardt fc Elvisl |2010| ) and II 
ISteinhardt fc Elvisl |2009al ) . quasars in the SDSS DR5 cat- 
alog i|Schneidenl200^ " are bounded at every redshift z < 4.1 
(well-defined at z < 2.0) by a sub-Eddington boundary 
(SEB) on the bright end and by a high-mass, low-luminosity 
boundary (HMLLB) on the faint end of the quasar sam- 
ple. These two boundaries restrict the observed luminosi- 
ties for quasars at fixed mass and redshift to a range 
spanning ~ 1 dex. These boundaries are defined by 95% 
dropoff in n umber density from the p eak number density at 
each mass ISteinhardt fc Elvisl [jpich . The location of this 
range evolves with both mass and redshift, lying within 
the boundaries given in Table [1] Because H/3-based and 
Mgll-based virial mas s estimators appear to be better than 
ClV-based estimators (jShen et al . 20081; iMarconT et al ]|2009l ; 
ISteinhardt fc Elvisl [2009^ 7 we will restrict ourselves these 
two estimators, and therefore to 0.2 < z < 2.0, when pro- 
ducing individual quasar tracks. 

For each of the quasars in the SDSS catalog, an ac- 
companying optical spectrum allows a determination of 
virial SMBH mass, bolometric luminosity, and redshift. Each 
SDSS quasar is one snapshot of the evolution of a particu- 
lar SMBH along some track in mass, luminosity, and time. 
The parameters describing that track describe the accretion 
history of the SMBH and reflect the feedback mechanism 
between the SMBH and its surrounding environment. Be- 
cause the final mass and path in mass, luminosity, and time 
taken by each model SMBH is sensitive to all four param- 
eters, successful model quasar tracks would correspond to 
a highly specific accretion history and therefore inform us 
about the nature of AGN-galaxy feedback. 



3 QUASAR ACCRETION HISTORIES: 
TRACKS IN THE MASS-LUMINOSITY 
PLANE 

We consider two models for the decline in Eddington ra- 
tio L/LEdd as a function of time: (1) an exponential decay 
L/L Edd = e ~ l - t ~ t °' >/k < ! and (2) a power-law decay L/L Ed d = 
(t — to)~ kp . The first model is inspired by arguments in 
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jHopk ins fc Hernauistll2009h . while the second uses a sim- 
ple alternative form to give an indication of the dependence 
on the assumed function. Both models are four-parameter 
fits, requiring a mass Mo and time io where the quasar de- 
parts from its Eddington luminosity, a decay slope k, and 
an accretion conversion factor ft = L/M — e/l, the ratio 
of the radiative efficiency e to the fraction of its lifetime I 
that the quasar is turned on (its duty cycle) . If the quasar is 
always turned on, the growth rate will be the instantaneous 
ft = L/M — e Our models make the simplifying assumption 
that the duty cycle and radiative efficiency are not time- or 
mass-dependent. Observations constraining quasar duty cy- 
cles, or a theoretical model constraining radiative efficiency, 
could allow this degeneracy to be broken. 

Existing observations and theory restrict the range of 
accretion conversion factors ft. We can estimate I from the 
observed fraction of galaxies hosting quasars fg. At low red- 
shift fg < 0.05, and it incre ases to approach fg ~ 0.2 
near z — 2 ("quasar evolution" . iMartini et all (|2009l i). Some 
SMBH might either not yet have reached a quasar state or 
have already turned off at the time fg is calculated. 

The observed fraction of quasars fg is the product of I 
and the fraction of galaxies going through a quasar phase at 
that redshift. Therefore, a measurement of fg — 0.2 places 
a lower bound of 0.2 on I at z — 2, but I might be con- 
siderably larger if relatively few galaxies contain a SMBH 
in a luminous accretion phase at that redshift. If only 20% 
of galaxies were capable of producing a quasar at that red- 
shift, quasars might be permanently on from the start of 
their quasar phase until turnoff. Theoretical models of accre- 
tion disks produce radiativ e efficiencies in the 8-2 0% range, 
with 10% perhaps favored |Yu fc Trem ainc 2002). Combin- 
ing these limits gives 0.08 < ft < 1, with the preferred values 
of I = 0.2 and e = 0.1 yielding a preferred value of ft = 0.5. 

For each set of model parameters, we construct the path 
in the M — L plane taken by one quasar obeying that model. 
Figure[T]shows the track taken by an individual quasar start- 
ing at 1O 8 ' 5 M0 and L^dd and declining in Eddington ratio 
as a power law L/LEdd = (t — to)~ 6 with ft = 0.5 (black, 
solid line). 20% increases in the initial mass (red), initial 
time (magenta), power- law slope (green), or k (blue) pro- 
duce different luminosities at higher mass, spanning ~ 0.6 
dex in total. The luminosity is most sensitive to changes in 
the power-law slope k and conversion factor ft. Changes in k 
and ft are nearly degenerate (steep values of k are preferred, 
as in §2J. 

We can also consider the mass and luminosity evolution 
as a function of time. In the resulting t — M and t — L planes 
(Figures[2]and[3|, changes in k and k (again, green and blue, 
respectively) now produce slightly different shapes, and the 
final SMBH mass is sensitive to each of these parameters. 
The exponential decay model exhibits sharper dependence 
upon track parameters in the t — M plane but weaker de- 
pendent in the M — L plane, as illustrated in Figures, [4][5l 
and [5] 



4 ALLOWED PARAMETER RANGES 

Using the two quasar track parametrizations from the previ- 
ous section, we consider how individual quasars can evolve in 
mass and luminosity. In this section, several methods are de- 
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Figure 1. Black: Evolution in the mass-luminosity plane of an 
individual quasar from 1O 8 5 M0, starting at the Eddington lumi- 
nosity and declining in Eddington ratio as L/L^dd = (* — *o) — 6 i 
or k p = 6. 20% increases in the initial mass (red), initial time (ma- 
genta) , power-law slope (green) , and conversion factor k (blue) are 
also indicated. The Eddington luminosity is shown as the black, 
dashed line. 



10.5 




6 8 
Time (Gyr) 

Figure 2. Black: Evolution in the time-mass plane of an individ- 
ual quasar from 1O 8 ' 5 M0, starting at the Eddington luminosity 
and declining in Eddington ratio as L/ L Edd = (t — to) - 6 - 20% in- 
creases in the initial mass (red), initial time (magenta), power-law 
slope (green), and conversion factor k (blue) are also indicated. 



scribed for constraining the evolution of individual quasars 
by requiring that their track lie within the observed M — L 
locus at every redshift, as in Figure [7] The simplest method 
is to use the boundaries from Table [1] to restrict allowed 
mass and luminosity as a function of time. We then proceed 
to consider further restrictions due to intrinsic quasar vari- 
ability, methods for dealing with digitization problems when 



4 Charles L. Steinhardt, Martin Elvis, and Mihail Amarie 




Time (Gyr) Time (Gyr) 



Figure 3. Black: Evolution in the time-luminosity plane of an 
individual quasar from 1O 8,5 M0, starting at the Eddington lumi- 
nosity and declining in Eddington ratio as L/L^dd = (t — *o) — 6 - 
20% increases in the initial mass (red), initial time (magenta), 
power-law slope (green), and conversion factor re (blue) are also 
indicated. Note that although the Eddington ratio is strongly de- 
creasing, the increasing black hole mass results in a luminosity 
that first increases, then decreases. 




Log M / M 

Figure 4. Black: Evolution in the mass-luminosity plane of 
an individual quasar from 1O 8 5 M0, starting at the Edding- 
ton luminosity and declining in Eddington ratio as L/L^dd 
■ or fc e = 1 Gyr. 20% increases in the initial mass 

(red), initial time (magenta), power-law slope k p (green), and 
conversion factor re (blue) are also indicated. The Eddington lu- 
minosity is shown as the black, dashed line. 



Figure 5. Black: Evolution in the time-mass plane of an individ- 
ual quasar from fO 8 ' 5 M0, starting at the Eddington luminosity 
and declining in Eddington ratio as L/L^dd = e_ t' - ' )/ 1 - , or 
k e = 1 Gyr. 20% increases in the initial mass (red), initial time 
(magenta) , exponential decay time k e (green) , and conversion fac- 
tor re (blue) are also indicated. 
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in the time-luminosity plane of 
8 - 5 M , starting at the Edding- 



Figure 6. Black: Evolution 
an individual quasar from IO'-'-mq 
ton luminosity and declining in Eddington ratio as L/L^dd = 
e — (*— *o)/ 1 -O i or A; e = 1 Gyr. 20% increases in the initial mass 
(red), initial time (magenta), exponential decay time k e (green), 
and conversion factor re (blue) are also indicated. Note that al- 
though the Eddington ratio is strongly decreasing, the increasing 
black hole mass results in a luminosity that first increases, then 
decreases. 
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Figure 7. The portions of two sample tracks (red) with Mo = 
lO 8 - 5 M ,t o = 3.5 Gyr, k £ = 1.0 and 0.5 Gyr, k = 0.65 and 0.35 
lying in each redshift bin (blue) are constrained to lie within the 
observed quasar mass-luminosity locus (purple boundaries). 



Figure 8. Allowed (green) combinations of slope and conver- 
sion factor for quasars with a power-law-decaying Eddington ratio 
starting at M = lO 8 5 M ,t o = 3.5 Gyr (z = 1.85). 



crossing from one redshift bin to the next in Table [TJ and 
methods for dealing with potential mismatches between the 
M — L plane as calculated using H/3-based and Mgll-based 
virial masses. 

Using the most straightforward of these methods, di- 
rectly applying the boundaries as in Table [T] consider a 
quasar at 1O 8 - 5 M and L E dd at 1.8 < z < 2.0, t = 3.5 
Gyr ( z = 1.85 using cosmological parameters from ISpergell 
(20071) ). Evolving this quasar forward, while constrained to 
lie within the Table [T] parameters, restricts the allowed slope 
k and conversion factor k to narrow bands for both power- 
law (Figure [5} and exponential (Figure decays. The al- 
lowed parameters form roughly a one-dimensional locus. As 
indicated by the test parameters in §[3l the turnoff time tj 
and final mass Mf are sensitive to each of these parame- 
ters. Requiring that the turnoff happen at a specific time or 
mass restricts the possibilities to a narrower set of possible 
parameters (e.g., 9.5 < log M/Mq < 9.6 in Figure [5]). 

Solutions with larger k, correspond to SMBH that gain 
mass more slowly at a given luminosity and/or are more 
often quiescent, resulting in slower growth and a smaller 
mass at turnoff. Solutions as k ~ s- rapidly produce high- 
mass SMBH that then turns off. These solutions therefore 
correspond to solutions in which growth is primarily non- 
luminous, and thus are disallowed by the requirement that 
most SM BH mass was acquired during a luminous quasa r 
phase (|Soltanlll982l ; lYu fc Tremaine]|2002l ; lElvis et alJl2002T ). 
and so are not plausible candidates for SMBH growth. 

4.1 Effects of Variability 

Individual quasars in SDSS are observed to undergo ~ 0.3 
dex variations in the optical on timescales of months to years 
l|de Vries et al.ll2005h . The luminosities on the model tracks 
represent the characteristic luminosity of the quasar aver- 
aged over longer timescales. However, the observed distri- 
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Figure 9. Allowed (green, magenta) combinations of e-folding 
timescale and conversion factor for quasars with an exponentially- 
decaying Eddington ratio starting at Mo = lO 8 ' 5 M ,tn = 3.5 
Gyr (z = 1.85). The magenta combinations correspond to quasars 
with a final mass of 9.5 < log M/Mq < 9.6. 



bution of quasars includes a dispersion due to this intrinsic 
short-term variability, the boundaries in Table [T] are sim- 
ilarly broadened, so the proper track must lie well within 
these boundaries. Individual quasars might spend a small 
fraction of their time across a boundary (5%, see § [2J ■ 

A restriction to lie 0.2 dex away from each boundary 
produces no non-trivial solutions at to = 3.5, Mo = 8.5. 
This might be because (1) the objects observed at 1O 8 ' 5 M 
are, due to intrinsic variability, above their characteristic 
luminosity, (2) the redshift bins are too wide, or (3) quasar 
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Figure 10. Allowed combinations of e-folding timescale and con- 
version factor for quasars with Eddington ratio declining both as 
a power law and an exponential starting at to = 3.5 Gyr and re- 
maining at least 0.2 dex from boundaries on the bright and faint 
ends of the quasar M — L distribution until turnoff . Initial masses 
are from 10 813 - 10 8 - 28 M Q . 

evolution does not lie on a simply-parametrized track of the 
sort we are fitting. If Mo is allowed to vary, keeping to = 3.5 
Gyr, solutions exist for the narrow range 8.13 < Mo < 8.28 
(Figure [TO). 

4.2 Common Tracks at all Mo, to 

One natural possibility is that every SMBH grows with the 
same value k and k but with different initial to and Mo- 
If so, this would be a strong indication that there is one, 
universal feedback process and that this simple model of 
quasar evolution has validity. Figures QT] and [12] show the 
allowed values of k and k for all power-law and exponential 
tracks, respectively, with any Mo from 10 5-10 , departing 
from Eddington at to = 3.5, 4.5, 5.5, and 7.5 Gyr. Here, we 
restrict solutions to lie within each boundary at all times, 
and in § 14. 31 will make these boundaries more restrictive due 
to quasar variability. 

We can also restrict the above solutions to lie at least 
0.2 dex from each boundary to account for scatter due to 
intrinsic variability in the observed quasar population. 

If all quasars follow a track with identical k and k, that 
track must be a valid solution in all panels of Figure [TTJ or 
1121 This would restrict possible combinations of k and k to 
the ranges shown in magenta in Figures [TTJ and 1121 There 
exist families of quasar tracks with the same k and k valid 
at all times to > 3.5 Gyr (z < 1.91). 

In addition to those shown in Figures [TT1 and [T2l there 
are also solutions involving quasars with small (Mo ~ 
10 6 M Q ) initial masses spending most of their lives at lu- 
minosities below the SDSS detection threshold and only be- 
coming massive and luminous enough to lie within the SDSS 
catalog with L <C Lsdd and near z = 0. Although these solu- 
tions do lie within the boundaries in Table [JJ such solutions 



Figure 11. Allowed combinations (green, magenta) of power-law 
slope k and conversion factor re for quasars with Eddington ra- 
tio declining as a power law starting at four different choices of 
to- The boundaries at each redshift are discontinuous, as given 
in Table [TJ Discontinuities in the boundaries on allowed combi- 
nation correspond to the quasar turnoff time moving from one 
redshift bin in Table [TJ to another. Magenta combinations of fc,re 
are allowed in all four panels. 
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Figure 12. Allowed combinations (green, magenta) of e-folding 
timescale k and conversion factor re for quasars with an 
exponentially-decaying Eddington ratio starting at four different 
choices of to- The boundaries at each redshift are given in Table 
[TJ Discontinuities in the boundaries on allowed combination cor- 
respond to the quasar turnoff time moving from one redshift bin 
in Table [TJ to another. Magenta combinations of fc,re are allowed 
in all four panels. 
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Figure 13. Allowed combinations of e-folding timescale and con- 
version factor for quasars with Eddington ratio declining as a 
power law starting at four different choices of to- The boundaries 
are interpolated from those in Table [TJ 
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Figure 14. Allowed combinations of e-folding timescale and con- 
version factor for quasars with an exponentially-decaying Edding- 
ton ratio starting at four different choices of to . Boundaries on the 
quasar M — L locus are interpolated from those in Table [TJ 



cannot describe the quasar distribution at z > 0, not do 
these solutions contribute to the Soltan argument (|Soltanl 
1982). We have only included solutions in both this and 
the following sections (Figures llll416|) that describe quasars 
which would lie in the SDSS catalog with Mbh = Mo and 
redshift corresponding to to. 

4.3 Other Interpolation Methods 

We can also restrict the above solutions to lie at least 0.2 
dex from each boundary to account for scatter due to in- 
trinsic variability in the observed quasar population. How- 
ever, there are no non-trivial solutions lying at least 0.2 dex 
away from each boundary with to = 4.5 Gyr. This may 
be because the boundaries change discontinuously when the 
quasar redshift moves across a bin boundary in Table [T] and 
some boundaries in this table match poorly. The greatest 
mismatches occur when virial masses transition from using 
Mgll to H/3 broad emission lines at z = 0.8. There are many 
ways one might correct for this binning problem. As here we 
are simply developing models in an effort to establish rough 
guidelines for what properties are required of quasar evolu- 
tion, we will just give the results of using a simple interpola- 
tion. In Figures [T3l and [T4l we take the boundaries in Table 
[T]to be correct at the midpoint (in time) of each bin, and in- 
terpolate in both slope and intercept to find the boundaries 
used at other redshifts. Using these interpolated boundaries, 
Figures [T3] and [TJ include parameters for all tracks of the 
sort included in Figures [TT1 and [T2l with 10 5 < M < 10 10 , 
departing from Eddington at the same four to. 

If quasar activity at all times is controlled by the identi- 
cal physical constraints, all quasars might lie on a family of 
tracks with identical shape (k and k) but different Mo and 
to- There are no quasar tracks with either a power-law or 
exponential decline in Eddington ratio with k and k lying in 
the allowed range in all panels of either Figure. Thus, fami- 



lies of tracks with the same k, k at all times are disallowed 
with this interpolation. 

4.4 Restriction to Mgll masses 

Finally, recent studies comparing H/3- and Mgll-based 
virial masses have proposed that a substantial correc- 
tion might be required to bring the two in to agreement 
i|Onken fc Kollmeierll2008l ; lRisaliti et al.ll2009l ). although the 
correction might not be n eeded for the brightest quasars 
|Steinhardt fc Elvis|[2009bh . The allowed parameters in the 
three panels with to < 6.0 Gyr of Figures ITTI and [T2l are pri- 
marily determined using Mgll masses, while the to = 7.5 Gyr 
(z — 0.68) panel is determined using primarily H/3 masses. 
Requiring that parameters appear only in all three Mgll 
panels produces the allowed parameters in Figures [Tol and, 
for interpolated boundary conditions, the power-law param- 
eters indicated in Fig. 1161 There are still no parameters for 
an exponential decay that appear in all three Mgll panels of 
Figure [12] (for the interpolated boundary conditions) . 



5 DISCUSSION 

We have shown that the quasar mass-luminosity plane can 
be used to restrict possible tracks for the evolution of in- 
dividual quasars. The allowed parameters presented in this 
paper correspond to simple models and analysis methods. 
In § [5] we discuss some of the properties that appear to 
be common to all allowed tracks. Additional tuning of the 
quasar track shape, e.g., requiring a different minimum dis- 
tance from boundaries on the quasar distribution, or the use 
of interpolation methods produces altered sets of allowed pa- 
rameters. These figures are intended to serve only as a rough 
guide to the types of changes that result from using different 
assumptions. 
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Figure 15. Quasar track parameters allowed for the three choices 
of tg < 6.0 Gyr (with Mgll-derived masses) in Figures 1 1 1 1 and [T2l 
The boundaries at each redshift are given in Table [l] 
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Figure 16. Quasar track parameters allowed for the three choices 
of to < 6.0 Gyr (with Mgll-derived masses) in Figures 1131 and 
1141 There are no tracks allowed at all three choices of to for an 
exponential decline in quasar Eddington ratio. The boundaries at 
each redshift are interpolated from those in Table [l] 

The production of these tracks was motivated by a con- 
sideration of the surprising synchronization of quasar lu- 
minosity evolution at fixed black hole mass. It was hoped 
that the seeming independence of quasar accretion rates on 
host galaxy parameters might allow quasar observations at 
different redshifts to be matched as lying along a common 
evolutionary track. While it may have seemed highly un- 
likely that the many complications provided by galactic dy- 
namics, merger history, and asynchronous quasar behavior 
were all truly negligible, surprisingly we find that simple, 



parametrized models indeed may provide a good description 
of the evolution of individual quasars. Further, it may even 
be possible that one family of similar tracks can describe the 
evolution of all Type 1 quasars. 

5.1 Tracks 

We have tested models with both an exponential and a 
power-law decline in quasar luminosity when constrained 
to pass through the observed quasar locus at all lower 
z. The power-law decline model has been proposed by 
iHopk fc Hernquistl |2009). The quasar tracks are sensi- 
tive to ~ 20% parameter changes (Figures[T]and[2}, showing 
that these tracks provide a sensitive probe of the details of 
quasar fueling and feedback. 

While the allowed ranges are sensitive to the details of 
how the boundaries of the quasar distribution are defined 
at each redshift, the allowed tracks have some characteristic 
general properties: 

(i) Typically, tracks are allowed over a restrictive range of 
rates of Eddington ratio decline and "accretion conversion 
factors" k. Different combinations of decline rate k and n 
result in different amounts of luminous mass growth from the 
start of the quasar phase until turnoff, as shown in Figure 

m 

(ii) All allowed tracks are characterized by a steep decline 
in Eddington ratio. For a power-law decline L/LEdd oc t~ k ° , 
even allowing n to be as large as 2 results in a track with 
fc e > 4, while placing k < 0.5 as expected would result in 
declines with k e > 6 or steeper. 

5.2 Rapid Decline 

We might have expected that only densities of gas and 
dust within the host galaxy density should be impor- 
tant, but the synchronizat ion of the quasar accretion rates 
jSteinhardt fc Elvid l2009ah implies that this might not be 
the case. In a matter-dominated universe (as existed at 
z = 2), the matter density falls as (1 + z)~ 3 oc t~ 2 . If quasar 
fuelling were linked to densities in the intergalactic medium, 
then, the falloff could not be t~ 6 or steeper. The rate of 
random two body collisions, e.g., galaxy collisions is propor- 
tional to the square of their number density in the universe, 
giving t~ 4 . A dependence of t~ e ' or steeper could arise from 
collisions involving at least three bodies. Alternatively, other 
astrophysic s might dictate a steeper decline, as proposed by 
iPena (| 20101 ) for "quenching" processes dominating galactic 
evolution. In the exponential case L/LEdd oc e~ h ", the de- 
cline must also be quite sharp, typically with an e-folding 
time of 1 Gyr or shorter. 

Both sets of allowed tracks present quasars as relatively 
short-lived objects, living only ~ 1 — 2 Gyr from their ar- 
rival above SDSS detection near the Eddington ratio until 
turnoff. The SDSS includes a large quasar population on 
0.2 < z < 4.1, a range spanning 10 Gyr. If the typical quasar 
can indeed only exist in the SDSS catalog for ~ 1 — 2 Gyr, 
then a typical galaxy can only contain a quasar a maximum 
of ~ 10 — 20% of that time, even if the quasar duty cy- 
cle includes no periods of quiescence. A mostly-quiescent 
duty cycle would further reduce the observed fraction of 
galaxies containing quasars. Since this observed fraction is 
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~ 10 — 20% near z — 2 and only slightly low er towards the 
edges of the SDSS quasar redshift range (cf. iMartini et al.l 
(2009)), we conclude that the typical quasar indeed has a 
duty cycle such that it is luminous closer to 100% than to 
10-20% of the time. For the duty cycle to fully reach 100%, a 
track with either a large (~ 25% or more) radiative efficiency 
or a steep decline (towards t~ 10 ) appears to be required. The 
implication is that SMBH evolution might not involve "flick- 
ering", i.e., interspersed periods of l uminous accretion and 
quies cence on a ~ 10 7 year timescale (|Hatziminaoglou et al.l 
1200 j ). 

Instead, a simpler picture suggests itself: (1) the SMBH 
is seeded; (2) the SMBH grows until it enters the SDSS 
catalog at the low-mass end and at the Eddington luminos- 
ity; (3) the SMBH accretes as a Type 1 quasar for 1 — 2 
Gyr while the Eddington ratio declines sharply; and (4) the 
SMBH permanently ceases its rapid, unobscured, luminous 
accretion, with various possible post-turnoff states including 
quiescence, Seyferts, and Type 2 quasars. This picture also 
requires a low value of k, corresponding to decline propor- 
tional to t -6 or steeper. 

The possibility of a short, 1 — 2 Gyr quasar life- 
time is particularly intriguing in light of two additional re- 
sults from the M — L plane. The characteristic luminos- 
ity for quasar accretion at fixed mass and redshift requires 
that accretion rates be sy nchronized to within 1 — 2 Gyr 
llSteinhardt fe Elvisll2009al ). Also, quasar turnoff is synchro- 
nized, depending upon the mass, to w i thin .75 — 3 Gyr for 
M B H > 10 9 M© i|Steinhardt fc Elvis! l2009ah . Perhaps the 
similarity of these three synchronization timescales could 
be explained by a synchronization in the times with which 
quasars turn and follow a common track in the M — L plane 
combined with short lifetimes. 

A short-lived Type 1 quasar phase might seem to vio- 
late the Soltan argument because the SMBH spends most of 
its time in another state. The Soltan argument shows that 
most of the total quasar mass in the universe was accreted 
luminously in Type 1 quasar states. However, the Soltan 
argument only places a limit on the last 2-3 e-foldings of 
mass growth, during which most of the mass is added. Prior 
to these last e-foldings, we do not know how much of the 
growth takes place through luminous accretion. Even these 
short tracks with decline between t~ 6 and t~ 10 grow the 
SMBH by 1-1.4 dex, i.e., 2-3 e-foldings, so such solutions 
are allowed. 



5.3 A Common Track 



As shown in Figures Figures [ill 121 [151 and 1 161 it is possible 
that one scaling law for feedback with universal parameters 
might be able to describe the evolution of all quasars at all 
initial masses and times. The existence of a characteristic 
luminosity at each combination of mass and redshift is in- 
sufficient to require such a uniformity among the evolution 
of individual quasars. However, (1) the synchronization be- 
tween quasars at fixed mass, (2) the narrow luminosity range 
at fixed mass and redshift, and (3) the sharp peak in number 
density at a sing le, characteristic luminosity at fixed mass 
and redshift (cf. ISteinhardt fc Elvis! l|2009ai )) make such a 
model intriguing. 



6 CONCLUSIONS 

In this paper, we have investigated tracks for SMBH ac- 
cretion histories and have shown that the quasar mass- 
luminosity plane constrains these models remarkably pre- 
cisely. The most intriguing result is that we can rule out 
models in which the SMBH accretion rate is proportional 
to the matter number density n in the universe. Even mod- 
els in which the accretion rate is proportional to n 2 are not 
allowed without a combination of higher-than-expected ra- 
diative efficiency and a SMBH that spends most of its time 
quiescent, rather than in a quasar state. 

This paper is an intermediate phenomenological step 
that has produced constraints that seem to be required of 
theoretical models for SMBH growth. The next step is to 
produce a physical model of fuelling and feedback leading 
to quasar tracks satisfying these constraints. 
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